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Abstract 

The electrons on the surface of a disordered multi-layer integer quantum Hall sys- 
tem constitute an unusual chiral metal with ballistic motion transverse to the field, and 
diffusive motion parallel to it. We present a non-perturbative analytic treatment of an 
' appropriate model, consisting of disordered chiral Fermions in two dimensions. A super- 

symmetric generating functional is set up for the correlation functions of this system. The 
C^h! strong disorder limit is mapped into a supersymmetric spin chain, with ferromagnetic ex- 

change coupling, reflecting the electron's chiral motion. The ferromagnetic ground state 
and the spin wave excitations, corresponding to the diffusion modes of the chiral metal, 
are found exactly. The parametric density of states correlator in the ergodic limit is com- 
puted from a Boltzmann- weighted sum over low energy spin states. The result is of a 
universal form and coincides with that for a Hermitian random matrix. 
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§ 1 Introduction 



Disorder has a profound effect on low dimensional electron transport, generically leading to lo- 
calization of all states in both one and two dimensions (2d) JTJ . In the quantum Hall regime with 
a strong magnetic field, however, the behavior is richer 0. Precisely at the transition between 
successive integer Hall plateaus, the 2d electron states are not localized, but quasi-extended 
with multi-fractal scaling characteristics 0. Within a Hall plateau, the Id edge states are ex- 
tended, due to their purely chiral nature[|, ||. Indeed, for a single edge mode, impurities simply 
.£h ! lead to an unimportant forward scattering phase shift. However, when multiple edge modes 
are present, such as at the boundaries of hierarchical fractional quantum Hall effect (FQHE) 
o3 ' states||, interchannel impurity scattering can be important 0. For example, backscattering 
between the two counter propagating modes of a v = 2/3 state is predicted to drive a transi- 
tion into a phase with charge and neutral sectors decoupled Multiple edge modes are also 
present in multi-layer systems exhibiting a 3d bulk QHE in perpendicular applied field. To- 
gether, these edge modes comprise a conducting two-dimensional subsystem, which has been 
the focus of recent attention|| 110]. Such systems can be realized by fabricating multi-layer 



GaAs heterostructuresfll|], but occur naturally in the Bechgaard salts, a class of quasi Id com- 



pounds which exhibit a cascade of field-induced spin-density-wave transitions between bulk 



QHE phases [12 



The chiral 2d metal "living" on the surface of such bulk QHE systems can be probed via 
transport experiments along the field direction (z-axis). For the integer quantum Hall effect 
(IQHE), the edge state in each layer is a free (chiral) Fermion, so that uniform inter-layer 
tunnelling leads to coherent z-axis motion. The low energy surface excitations then comprise 
a 2d Fermi surface, which is equivalent to "half of a conventional 2d open Fermi surface. A 



standard diagrammatic approach[10] shows that impurity scattering leads to diffusive motion 
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along the z-axis, with ballistic motion perpendicular to the field. Since backscattering is not 
possible in the ballistic direction, conventional localization effects are not expected. Indeed, 



a perturbative expansion about the diffusive metal shows an absence of any localization [10 



Moreover, random- walk arguments and numerics on an appropriate network model suggest 
that this conclusion is valid generally. 

In this paper we describe a non-perturbative analytic treatment of the IQHE surface sheath 
with impurity scattering. We show that an appropriate generating functional for the 2d electron 
correlations can be mapped onto a Id quantum spin chain. The replica trick is avoided by 
employing Fermionic and Bosonic partners (Bosonic "ghosts" ) from the outset to perform the 
desired quenched average. The resulting spin chain is thus supersymmetric, involving both 
Bosonic and Fermionic spins. Since the edge modes all move in the same direction, the super 



spin chain has a ferromagnetic exchange interaction [y|]. We show that the ground state of this 



super spin chain consists of all spins "aligned", analogous to the ground state of an ordinary 
ferromagnet. Moreover, there is a class of one-magnon excitations which are exact eigenstates of 
the super spin chain, leading to a sharp one-magnon pole in the spin-spin correlation function. 
This pole corresponds to the diffuson pole, describing the diffusive z-axis motion of electrons 
in the chiral 2d system. The exact pole confirms the complete absence of all localization 
corrections for this 2d chiral metallic system. 

The use of supersymmetry has a long history in the theory of disordered electron transport. 
Efetov[|14| reformulated earlier replica field theories |U| to obtain a supersymmetric non-linear 
sigma model (NLoM) description of disordered electron motion. This approach enabled the 
computation of several new properties, most notably electron localization in Id and random 
matrix spectral correlations in Od. To describe plateau transitions in the IQHE, a topological 
term is needed in the cr-model, as originally argued by Pruisken et. al.|16|] and incorporated into 
the supersymmetric approach by Weidenmuller fP7|| . More recently, one of us|18j has shown that 
the 2d supersymmetric cr-model with topological term can be mapped into a Id super spin chain 
Hamiltonian, with antiferromagnetic exchange interaction. This spin chain was, in turn, shown 
to be equivalent to the Chalker-Coddington network modellT^, EUf. The antiferromagnetic 



exchange corresponds to counter propagating edge modes of the network model, in contrast to 
the ferromagnetic exchange for the chiral 2d metal studied here. Unfortunately, the AFM super 
spin chain model is much more complicated than the ferromagnetic model analysed below, and 
has so far eluded analysis. 

In addition to obtaining the exact ground state and one-magnon excitations, we employ 
the ferromagnetic super spin chain to extract spectral correlations of the chiral metal. A finite 
length spin chain at non-zero temperatures corresponds to a 2d chiral metal with both dimen- 
sions finite. As with other zero- dimensional random electron systems, we expect interesting 
universal spectral correlations ||2l|1. Working directly with the Id spin Hamiltonian, we derive 
expressions for the two-point spectral correlation function and the parametric correlations. Not 
surprisingly, the results coincide with the universal random matrix theory results for the uni- 
tary ensemble, obtained previously from other methods fl22|, |23| . Our approach is notable for 
its (relative) simplicity, involving a simple sum over a large (super-)spin multiplet of quantum 
states. 

This paper is organized as follows. After introducing the model in Section II, we construct 
the supersymmetric generating functional in Section III, and show that in the limit of strong 
disorder it can be mapped into a super spin chain, with ferromagnetic exchange. In Section 



2 



IV we consider the electron's motion in the thermodynamic limit. From the exact ground 
state and one-magnon spectrum of the spin chain, the diffuson (or density-density correlator) 
can be computed exactly, revealing diffusive motion parallel to the field, and ballistic motion 
perpendicular to it. In Section V we consider mesoscopic effects for a sample with finite 
dimensions. Due to the system's chirality and anisotropy there are three mesoscopic regimes. In 
the "zero-dimensional" ergodic regime, we compute the two-point spectral correlation function 
and the parametric correlations, by summing over a zero-energy multiplet of super spin states. 
Section VI is devoted to a summary and brief conclusions. 



2 Model 




Figure 1: Geometry of a 3d quantum Hall sample, z-axis transport is included via the tunneling 
amplitude t, and impurity scattering by the random potential V. The circumference is denoted 
(3 = L x , and the height in units of the layer spacing is N. 

The system of interest (see Fig. 1) is a collection of N edge states, each described by a chiral 
Fermion, which are coupled together by an inter-edge hopping strength, t: 

N rfi 

n = J2 dx Wn id xi>n ~ *(Wn+l + h.C.)], (1) 

where (3 = L x is the circumference of the surface sheath (and the sample). Our choice of units 
is such that the edge velocity v — 1. The full Hamiltonian including random scattering is 
H — Hq + Hi, where 

Wi = E / dxV n (x)^ n . (2) 

n J 

For simplicity the random potential is assumed to have a Gaussian distribution with zero mean 
and 

[V n (x)V n/ (x')] cns = 2u5 nn/ 5(x - x'), (3) 
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where the square brackets denote an ensemble average over disorder realizations. 



3 Generating functional 

As usual, a generating functional can be introduced from which the electron Green's functions 
can be extracted Of interest are products of retarded and advanced Green's functions, 
denoted G±, where 

G±(n, n 1 ] x, x'\ ±u>/2) 
= (nx\ [H±i(ji + iuj/2)Y 1 \n'x'). (4) 

Here rj is a positive infinitesimal, and u/2 is the electron energy. The appropriate generating 
functional can be written as an integral over a spinor of Grassmann fields, ip na , where a = ± for 
retarded (advanced). To make the analogy with spin, we will sometimes use instead the notation 
a =|, j. To avoid the use of the replica method, we adopt the supersymmetric approach, and 
consider an additional integration over a spinor of complex fields, <p na . The supersymmetric 
generating functional, Z, introduced below, is then normalized, with Z = 1 for every realization 
of the disorder potential. An average over an ensemble of disorder realizations can then be safely 
performed. 

To insure convergence of the functional integrals over <p, the contribution to Z in the retarded 
sector is defined by exponentiating i(Ti. — u>/2 + irj), whereas in the advanced sector, —i(TC + 
uj/2 — irj). The appropriate supersymmetric generating functional then takes the form 

Z = J D*/jD^D(j)D(f)*exp(-S) } (5) 
with S = dx^2C n . (6) 

J n 

The Lagrangian is expressed as a sum of Boson and Fermion contributions: 

£ n = L(0*,0) + L(^), (7) 

with 

L^A) = <l>la z {d x -iV n )(l> n 

+ ^(0>A+i + cc)+#;0 n , (8) 

where r) = rj + iuj/2. In the future, we will drop the tilde, remembering to let rj —>■ rj + iuj/2 when 
necessary in computing quantities at non-zero frequency. By construction, Z = 1, independent 
of the random potential. 

The single particle Green's functions can be expressed as averages of either the complex or 
Grassmann fields. Specifically, 

G a (n,n';x,x') = —ai < (f) na (x)(f)n> Q (x') > (9) 
= -ai < ^ na {x)i) n , a {x') >, (10) 
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where a = ± and the brackets denote an average taken with weight exp(— S). The density of 
states (DOS) follows from the imaginary part of the diagonal element, 

N (3 

p = 7T- ( G - ~ G+)(n,n;x,x). (11) 
„=i 2m 

Note that an integral over x is unnecessary, since the DOS is independent of x. This follows 
from the chiral conservation law, 

d x \$\ 2 = V n - J, (12) 

where the current is J(x, n) = 2tlm[$*(a;, ra)$(x, n+1)]. Eq. ( p72|) is valid for every eigenfunction 
$ of the Hamiltonian. To resolve the ambiguity which arises because the DOS involves G at 
equal points, x = x', we take the symmetrized form, 

Anp = (3j2J2<(f) na (x±e) ( j)* na (x)> (13) 

n,a ± 

= PT,Y, < 'Pna(x±e)^ na (x) >, (14) 

n,a ± 

with e a positive infinitesimal. 



3.1 Coherent State 



Following D.H. Lee|K]], we reinterpret the spatial coordinate x as an imaginary time coordinate 



r, with 5* the Euclidean action for a Id quantum system. The Hamiltonian for this Id quantum 
system, denoted H, acts as a transfer matrix in the x direction. In order to extract H, it is 
convenient to recast the action S into the canonical form for a coherent state path integral 
of a Id system of Bosons and Fermions. Up to a sign in the advanced sector, the linear x- 
derivative term in (8) is already in the appropriate form. To "correct" this sign, we transform 
the advanced complex fields as 

<f>ni «-> Ki* ( 15 ) 

leaving the retarded fields, </> n f, unchanged. The first term in the Lagrangian then takes the 
canonical form: <fi* n d T (j) n . In the Fermion sector we similarly transform the Grassmann fields: 

^nl -> -VV> ( 16 ) 
VVil ^H, (17) 

again leaving the retarded Grassmann fields unchanged. The transformed Lagrangian can be 
written C n = L + L t + L v + L v , with 

L = (f)* n d T (f) n + iJ n d T tp n (18) 
L v = -iV n {r){4)* n a z (j) n + ^ n cr>J, (20) 



in canonical form. The term 
and the random potential, 
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are unchanged, whereas the tunnelling terms become 



L t = it(A n+hn + A n>n+1 ) (21) 

with 

A m , n = CAt - 0; i mi + (0 - V)- (22) 
3.2 Id Hubbard model 

We can now perform an average over the ensemble of random potentials. This gives 

L u = u (f n o z $ n + tf n <T g ip n ) 2 ■ (23) 

The full action S is now equivalent to the Euclidian action for a Id quantum system, and the 
Id quantum Hamiltonian (transfer matrix) can be readily extracted. One simply replaces the 
Grassmann fields ip n by Fermion operators, /„, and the complex fields, <f) n , by Bose operators, 
b n . The resulting Id Hamiltonian takes the form: H — H t + H u + H v with 

H v = vT,( b lK + flfn), (24) 

n 

H u = uJ2(bla z b n + fla z f n ) 2 , (25) 

n 

(26) 

n 

where 

A n+1 , n = b[ +1 ^ - 6^6 n+li + (b - /), (27) 
and the Fermion and Boson operators satisfy canonical commutation relations 

lha, blp) = $nn>5 a p, (28) 
[L. flp] = Snn'Sap- (29) 

Here, and in the remainder of the paper, [0i,0 2 ] denotes the graded or super- commutator, 
defined as 

[0 ± , 2 ] = 1 2 - {-l)\°^0 2 O u (30) 

where \0\ = if O is a Bosonic operator, while \0\ = 1 if it is a Fermionic one. This Id 
Hamiltonian describes spinful Bosons and Fermions hopping on a Id lattice, with hopping 
strength t, interacting via a quartic term with strength u. Large interaction u corresponds 
to strong disorder. For large disorder u, this term constrains the on-site Bose and Fermion 
densities to satisfy 

&W + / T t / T = &fo + /|/ i , (31) 

on each site n. There is a large energy cost (u) for configurations involving different numbers 
of particles in the retarded and advanced sectors, whereas retarded and advanced particles can 
move together. This reflects the phase cancellation between the retarded and advanced Green's 
functions. The positive infinitesimal restricts the total number of particles. 
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The similarity between the full Hamiltonian and more familiar interacting Id models such 
as the Hubbard model |24j] can be revealed by introducing a new set of operators. To this end, 
we define 

£? T = 6 T , B i = &{, (32) 

and 

*T = /T> F l = fl (33) 
The F operators are bona fide Fermion operators, satisfying 

[^na,4/?]=<W<W> (34) 

but the operator does not satisfy the canonical Boson commutator, but rather, = 
— 1. To restore the canonical form we define 

B l = -B\, (35) 

and also B^ — Bt, which now satisfy 

[B na , B n ip] = 5 nn iS a /3. (36) 

As defined, B n and B n satisfy canonical Bose commutation relations, but it must be kept in 
mind that B n is not the adjoint of B n . 

In terms of these new operators, the Id Hamiltonian takes a form which closely resembles 
the Hubbard model, with spin-independent hopping and interaction: 

H t + H u = it ^2(A n+ltU + A njn+ x) 

n 

+uJ2(B n B n + FlF n ) 2 , (37) 

n 

where 

A n +i >n = B n+ iB n + F n+1 F n . (38) 

The positive infinitesimal term now acts like a small magnetic field, coupling to the ^-component 
of the "spin": _ 

H ri = V (B n a z B n + F^ z F n ). (39) 

As with the Hubbard model, one anticipates that the interaction term will lead to a "charge 
gap", strongly suppressing fluctuations in B n B n + F^F n . Of interest are the remaining gapless 
"spin" excitations. For the conventional Fermionic Hubbard model, these spin excitations can 
be revealed by employing a transformation, to map the Hubbard model into a Heisenberg spin 
chain A similar transformation is desirable in this Id super symmetric model. To this end, 
we will return to the path integral representation, and perform a gauge transformation before 
averaging over disorder. This will allow us to obtain directly a (supersymmetric) spin chain 
model. 
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3.3 Gauge transformation 

Consider the gauge transformation 

4> n -> exp(ia z ( dT'V n (r'))(j) n , (40) 

J — oo 

and an identical transformation for ip n . Due to the linear derivatives in Lq, this eliminates 
the on-site random potential term in Lv, at the expense of introducing randomness into the 
hopping term. The full Lagrangian becomes C n = Lq + L t + L v , with Lq and L v as given in 
Eq. (0) and Eq. (0), and 

L t = i(t n A n+ltn + t* n A Hin+1 ). (41) 
The tunnelling amplitudes, t n (r), are both random and complex, given by 

t n (r) = texp (i J^dr'[V n (r') - K+i(t')]) ■ (42) 

They satisfy 

[tn(T)t:,(T% ns = 5 nn/ t 2 e- 2u ^- T \ (43) 

In the following we will focus on the limit of strong disorder (large u), and approximate the 
above exponential with a delta function: 

[t n (r)t* nl (r% ns -> D8 nn ,5{r - r'), (44) 

where D = t 2 /u. The quantity D will play the role of the exchange interaction in the super- 
symmetric spin chain. The above large u limit is then analogous to passing from a Hubbard 
model to a spin chain. 

To extract the Id (spin chain) Hamiltonian, we again perform an ensemble average over 
disorder. Since the above delta function is really a short-range symmetric function of its 
argument, we will arrive at a symmetrized form. On passing again to Bose and Fermi operators, 
the full Id Hamiltonian is given by H — Hp + H v where H v was defined in Eq. (p9|) , and 



* n 

with A n+ i tn as before, in Eq. (|38|) . 

Before recasting the Hamiltonian if/) into the form of a spin chain, it is instructive to re- 
express the DOS in terms of the Fermion and Boson operators. Recalling Eq. ( |13D we obtain 
for the density of states, p = J2n Pn with 

4vrp n / (3 = (b^b n1 + b^b^ + b ] n[ b ni + 6^6^) , (46) 

where (...) is a "thermodynamic" expectation value with Hamiltonian H, see Sec. ^3] for the 
details. A similar expression exists in terms of the Fermion operators b —> f. Equivalently we 
have 

2np n /P = l + (blb n ) = l-(flf n ), (47) 

which involves the on-site Boson and Fermion densities. These expressions can also be written 

as 

2n Pn /f3 = (B n a z B n ) = - (F^a z F n ) . (48) 

As we shall now show, these correlators can be naturally interpreted as the z-component of an 
appropriately defined spin operator. In terms of these spin operators, the Hamiltonian Hd will 
take the form of a ferromagnetic (super-)spin chain. 
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3.4 (Super-)spin chain 

To cast Hjj in spin chain form, it is useful to introduce generalized spin currents. This can be 
done in a supersymmetric manner, by defining a four-component superfield, 

* n =(F n ,B n ), * n =(F*,B n ). (49) 

When necessary, we will use latin indices (a, &,•••) to denote the Fermion/Boson label, i.e. 
a — F, B <-> 1, 2. The \1/ obey the mixed statistics relation 

^maa^nbp = M ab ^ nbf3 ^ maa + S mn 5 ab 5 a /3, (50) 

where M ab = a x ah — a z ah . We may then define a four- vector super spin matrix via 

J& = ^al^ b = ^aal^ b p, (51) 

where 7 = (1, cr)/2. In matrix form, this is 

7" - ( F ^" F Fh " B ) (52) 

J ~ { B-fF B-fB J ■ 1 ' 

The global current operators 

^tot = ^ $n (53) 

n 

generate symmetries of H D , as can be easily seen by checking [j7"rp OT , y4 n+l n ] = using Eq. (|5Q|) . 
Thus 

[J&yr,H D ] = 0. (54) 

It is now straightforward to obtain the spin chain representation. Noting that the hopping 
operator A n+ i <n = \l> n+1 \l> n , we see that substitution in Eq. (|4"5]) yields 

n 

n 

The identities H nao $nbp = 2J^ a and Tr( 7 'V) = ^ u /2 then give 

H D = 2DJ2 Tr [j* - JfaJ^] • (55) 

n 

The minus sign on the spin-spin interaction indicates that this is a ferromagnetic (super-) spin 
chain.Q The 77 term indeed acts as an infinitesimal field 

H V = 2 V J2^J:. (56) 

n 

The term J2n Trj7^ = JTtot vanishes on the Bose- Fermi vacuum, and commutes with Hp and 
H v . It is therefore zero on all low energy states coupled by the Hamiltonian Hp + H v , and will 
be dropped from all future expressions. 

1 Note that the presence of a z is equivalent in the second term to the conventional "supertrace" used in 
supersymmetric notation. 
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In order to perform calculations, it is useful to write out the interactions in a more explicit 
(but not manifestly supersymmetric) form. Dividing the terms into three sectors, Hp = Hp + 
Hp + Hx, one finds 

Hf = -EWViAi + ^^Sj, (57) 

Z n 

= "tt ^ (B n B n B n+1 B n+ i + 4Sf ■ , (58) 



2 „ 
D 



n+1 



-B n F n Fl +1 B n+1 - A~S X n ■ S* +l ) , (59) 



where 



S n = F^F n , S* = B n ^B n , (60) 

&n - F ] n -B n , S n = B n —F n . (61) 

The Fermionic Hamiltonian, Hp, is precisely that of an ordinary spin-1/2 ferromagnetic spin 
chain ||24| . In particular, the Fermionic currents obey the SU(2) algebra, 

[SZ,S&=it*»SZ6 nn . (62) 

Since Hp is invariant under global spin rotations, we expect that the total spin operator 

Stot =J2 S n (63) 



n 



commutes with Hp. In fact, Eq. fl54j ) implies 

[5 TO t ,H f ] = 0= [5 TO t , H D ] , (64) 

i.e. the full Hp, is SU(2) invariant. H v of course breaks the symmetry explicitly. 

Although the Bosonic sector described by Hp appears also to be SU(2) symmetric, it is 
not so in the usual sense. Because of the unusual definition of conjugation (B ^ £?'), the 
currents are not Hermitian. This non-hermiticity means that, even though there exists a 
set of currents forming an SU(2) Lie algebra and commuting with Hp, these cannot be used to 
generate unitary transformations expm • S B . A Hermitian set can be defined by 

J x n = iSS u = {b^ + b n ^b nl )/2 i (65) 

jy = -iS** = {bl^ nl -b n ^b ni )/2i. (66) 

J z n = 5? = (6i T 6 BT + & B A)/2. (67) 
These obey instead the SU(1,1) algebra [p5|1 

[J* n ,Jl] = -Ul (68) 
[jy,J*\ = U:, (69) 

[J Z n,J X n] = ^1 (70) 
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where the minus sign in the first relation distinguishes SU(1,1) from SU(2). In terms of these 
operators, 

H B = 2L>^(i(6taA-l)(fcUi^A + i-l) 

n 

_l_ T z T z — T x T x — T y T y \ f71 

i u n u n+l u n u n+l < J n J n+lJ- \ IL ) 

The Lie algebra of SU(1,1) is a well-studied semisimple Lie algebra[|5|]. Quite generally, such 
Lie algebras have a so-called quadratic Casimir operator A, which commutes with all elements 
of the algebra, and is the analog of S ■ S for SU(2). For SU(1,1), 

a = (rf - (rf - (jyf. (72) 

It is straightforward to verify that [A, J] = 0, and that therefore the global generator 

j TOT = j n (73) 

n 

commutes with Hb- As before, in fact, 

[J TOT ,H D ] = 0, (74) 
as can be verified by explicit calculation. 

3.5 Boundary Conditions 

In the above development, we have implicitly assumed that the coordinate x, which runs along 
the chiral modes, is infinite. Since x has been replaced by imaginary time r, this corresponds to 
zero "temperature", (3 — oo, for the Id (super-)spin chain. Consideration of mesoscopic effects 
requires a finite system, with both N and (3 finite. Since the chiral modes cannot end, when 
j3 is finite it is necessary to consider periodic boundary conditions in the x (or r) direction, as 
depicted schematically in Fig. 1. The boundary conditions in the direction (n) transverse to the 
chiral currents, can be taken either open, corresponding to a cylinder, or periodic, corresponding 
to a torus. 

When the system is periodic in the x-direction, care is needed in defining the boundary con- 
ditions on the complex and Grassmann fields. Specifically, the generating functional, Z, will 
only equal one if the determinant from the Grassmann integration precisely cancels the inverse 
determinant from the integration over the complex fields. This requires that the boundary 
conditions on ip and be identical: 0(0) = (f>((3) and -0(0) = i/j((3). However, in the standard 
imaginary time Grassmann path integral representation for Fermions, the Grassmann fields 
satisfy anti-periodic boundary conditions, ij){0) = —ip(0). This discrepancy prevents us from 
making the transition from the path integral to the quantum Hamiltonian by simply replac- 
ing Grassmann fields by Fermionic operators. To remedy this, we perform a global gauge 
transformation, 

i(x)^e ! "/> n (x), (75) 

in the Lagrangian Eq. (||), which transforms the boundary conditions for the Grassmann fields 
from periodic to anti-periodic, and puts the Grassmann path integral in the standard form we 
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desire. Under this transformation the Lagrangian picks up an additional term of the form 



Uc = l -^^ncr z ^n- (76) 



Since the transformation Eqs. (1^-17) leaves this form unchanged, this corresponds to an addi- 
tional term in the Id Hamiltonian: 

H bc =%Y,fl°Jn- (77) 

P n 

Equivalently, what we need to do when computing thermodynamic averages with the Boltzmann 
weight exp(—/3H), is to take instead of the usual trace the supertra.ee: 

(O) = STre^O = Tr(-l) Np Oe^ H , (78) 

counting the contributions from states containing an odd number of Fermions, Np, with a 
minus sign. This latter procedure is the one adopted in the following. Thus we do not modify 
the Hamiltonian, but will remember to insert a minus sign for all states with odd Fermion 
number. 

4 Diffuson 

In this section, we calculate the chiral diffuson in the thermodynamic (/?, N — > oo) limit. It is 
defined by 

K n (x;r]) = \G + (n + n', n'\ x + x', x'\ 0)| 2 

L J ens 

= [G + (n,0;x,0;0)G_(0,n;0,x;0)] ens . (79) 

It may be calculated in several ways, depending upon which fields (0 or ip) are used to generate 
the Green's functions G±. In particular, 

K n (r) = (T T S-(r)S^0)) (80) 

= <T t J-(t)J+(0)> (81) 

= -(%A-{r)A+m (82) 

= (TrA-(T)A+(0)), (83) 

where 

S n = Sn^iSni A n =S* x ±iS n y , (84) 
These expressions have a simple interpretation when expressed in terms of the b and / operators, 

&n = fn]fn[i = fn^n\l (86) 

■£ = 1>IA> A+ n= h Ufll (87) 
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The operators create pairs (for the Fermions, these are singlets) of particles at a particular 
site. Eqs. (|S0|-|8~3|) indicate that the diffuson is obtained by creating such a pair at the origin 



and propagating it to site n in time r. The motion of a pair describes the long-range coherent 
propagation, in the chiral disordered metal, of one retarded and one advanced particle. 



As discussed in Sec. p3j the angular brackets in Eqs. (p0|-|83|) indicate an average calculated 
by a trace with the weight exp(—f3H). In the thermodynamic limit, f3 — > oo, with 77 non- 
zero, only the zero energy (ground) states survive. Since the Hamiltonian is a sum of positive 
semi-definite operators, each of these operators must annihilate any prospective ground state. 
Satisfying this for the two terms in H v in Eq. (|24|) therefore implies that the only state which 
contributes to the average is the b — f vacuum |0), defined by 

b n \0) = f n \0) = 0. (88) 

In the Fermionic sector, this is nothing other than the fully polarized spin-down ferromagnetic 
state, so that 

SW|0) = -f |0>. (89) 
It follows that (as can be easily verified) 

S TOT \0) = 0, (90) 
i.e. |0) is the "lowest weight" state of this SU(2) (sub-) representation. Similarly, one has 

J TOT |0) = 0, (91) 
4ot|0) = y|0>, (92) 

so this is also the lowest weight state of the SU(1,1) (sub-)representation, and it is the SU(1,1) 
analog of the fully polarized state. As in the case of the usual ferromagnet, this is the exact 
ground state. 

Now consider the diffuson. From Eq. (BUI), we have 



K n (r) = <0|S-e-^S+|0> (93) 

= I ^(0|^e-^|0>e^, (94) 

where the "single magnon" operator 

St = E S ne ±tkn (95) 



simply creates a superposition of local spin-flips with wavevector k. To evaluate Eq. (|94]), we 
should study the state 

\k) = St\0). (96) 

Because no b Bosons have been added to the system, this state is still annihilated by H B and 
Hx- The problem then effectively reduces to a purely Fermionic one, i.e. the usual spin-1/2 
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SU(2) spin chain. As in that case, the state \k) is an exact eigenstatel By direct computation, 
one finds 

H\k)=E k \k), (97) 

where the dispersion relation is 

E k = 2.0(1 - cos k) + 2r] ^ 2r] + Dk 2 , (98) 

for k <C 1. This is just the usual ferromagnetic spin- wave dispersion. It is thus straightforward 
to evaluate Eq. (0), to obtain 

KJt -+x) = 9(x) / —e-^W- (99) 
jo 2% 

Performing the Fourier transform on x gives the final (exact) result 

K[k x , k z ) 



2t] + ik x + 2L>(1 - cos k z ) 

~ 2rj + ik x + Dkr (1 ° 0) 

We thus see that the ferromagnetic dispersion E k ~ Dk 2 + 2r) corresponds to the expected form 
of the anisotropic diffuson for the chiral metal ]10 . 



It is interesting to note that this result must also be obtained from the other formulations, 
Eqs. ( |3T| - |S3| ). In fact, one finds that there are three other sets of exact single super spin flip 
excitations, created by Fourier transforms of the operators J^~, A* , and in Eqs. 
The reader may amuse him or herself by carrying through this calculation explicitly. 



5 Mesoscopic Regime 

The above results describe the propagation of density fluctuations in the thermodynamic limit. 
We now consider the level statistics of a mesoscopic sample. A natural probe of these statistics 
is the density of states correlation function |23|j. 

CV; (3, N) = [p(E + uj/2)p(E - uj/2)] cns - f. (101) 

Of interest is the structure of this correlation function when the frequency u is on the scale of 
the mean level spacing, 

(102) 

where (3v is the system size in the x-direction, and v is the edge mode velocity. For isotropic 
metallic samples, C(u) is a universal function of the ratio lu/A, described by random matrix 
theory [P^, In the present case, the sample is both chiral and anisotropic, so care is needed to 
specify the random matrix theory regime. To this end it is useful to consider several important 
length scales. 
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5.1 Length scales 



For finite (3, the electron takes a finite time, t# ~ /3, to ballistically traverse the sample in the 
^-direction. Equating this with the time taken to diffuse across the system in the transverse 
direction, ~ N 2 /D, gives a (transverse) crossover length 



L = sJD(3, (103) 

measured in units of the z-axis lattice constant. This length scale has a simple interpretation 
in terms of the spin chain. For a chain with N sites, the energy of the lowest lying spin wave 
excitation, which has wavevector k m i n = 2n/N, is E min ~ D(k^ nin ) ~ r^ 1 . This energy equals 
the spin chain "temperature", E min ~ T — 1/(3, for a spin chain of length L . The boundary 
N = L thus demarkates the two following regimes: For a short spin chain with N < L , 
finite wavevector spin-wave excitations will not be appreciably excited, and the system will 
be effectively "zero- dimensional" . The opposite limit, N > L , which we refer to as a "one- 
dimensional diffusive" regime, corresponds to the "classical limit" of the spin chain. The 
transverse electron motion is diffusive. 

But there is another important length scale |§. A Id spin chain at finite temperature will 
have a finite spin-correlation length, 

£~D(3, (104) 

with D the spin-exchange constant. This length corresponds to the electron localization length 
along the z-axis. In the (3 — > oo limit of interest, one has £ ^> L . The above "Id diffusive 
regime" requires £ > iV > Lq. When N > £ the system enters into a third regime - a Id 
localized regime. 

Of the three regimes, the "zero-dimensional" is the simplest and the one we focus on. This 
regime corresponds to random matrix theory, as verified below. Since we consider frequencies 
on the scale of A, the Heisenberg time oj~ 1 ~ A -1 ~ N/3. In this regime, the ergodic condition 
that uj~ 1 exceeds both r B ~ (3 and t d requires simply N ^> 1. In the "Id diffusive regime", 
corrections to random matrix theory are anticipated, with wavefunctions possibly exhibiting 
multi-fractal scaling behavior, as in other low- dimensional diffusive systems @, [2(J. 

To evaluate the spectral correlator C(u) at finite (3 requires performing a full thermody- 
namic trace over states of the supersymmetric spin chain - a complicated task. Fortunately, 
in the "zero-dimensional" regime, all finite wavevector spin-wave excitations are exponentially 
suppressed (by e~ l3Emin = e~( 27rL °/ Ar ) ), and can be ignored. The trace can thus be restricted to 
a subset of states which are annihilated by the Hamiltonian Hp. This restricted trace is the 



operator analog of the zero-mode integration carried out on the conventional NLcrM[I3. 



5.2 Zero energy multiplet 

Consider then the zero energy multiplet of states that are annihilated by Hp, i.e. the set of 
ground states in the absence of the symmetry-breaking field. This set has a large degeneracy, 
because any particular vacuum (e.g. |0)) may be "rotated" using the symmetry operators 
J7xot t° obtain another zero energy state. These operators and the resulting multiplet form 
a representation of the U(2|l, 1) (super) symmetry group of H D , which contains SU(2) and 
SU(1,1) sub-representations generated using the operators Stot and Jtot- 
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To describe the E = multiplet, it is convenient to use the b Boson and / Fermion operators, 
for which |0) is the vacuum. Of the 16 elements of J~, 4 may be chosen diagonal. A natural 
choice is J® x , J X1) J® 2) J221 or > equivalently, the set of number operators, 



n 



fa 



^2,fnafnai n ba — E Hia^na- (105) 



All elements of J commute with the total b — f spin operator, which is 

Tr ^TOT = \ ( n n - n fl + n ^T - n bl) = 0, (106) 

where the final equality follows from evaluation in the |0) state. More fundamentally, this 
results from the Hubbard constraint Eq. ([31]). 

The remaining 12 elements of J act as raising and lowering operators. From Eq. (|54|) , we 
may act with any string of these operators upon |0) to produce another zero energy state. The 
12 operators, however, are strongly redundant, and in fact only 4 suffice to generate the full 
E = multiplet. 

The first two of these are simply the raising operators in the SU(2) and SU(1,1) subalgebras, 
which create delocalized "spinless" pairs of Fermions and Bosons, respectively, i.e. 

SioT = £/Uti> ( 107 ) 

n 

4ot = E & !A- ( 108 ) 

n 

These are just the k = magnon operators of the previous section. In addition, we must 
consider the Fermionic generators, 

^ot = E/!A> ( 109 ) 

n 
n 

These create Fermionic k = pairs. To simplify the notation we now drop the bulky subscript 
TOT. 

To see that S + , J + , A + , and A + , are sufficient to generate the full E = multiplet, note first 
that the remaining 8 off-diagonal generators all annihilate the vacuum (because they contain 
either b or /). Therefore, in any string of generators, they may be commuted or anticommuted 
to the right until they annihilate the vacuum, and it is straightforward to see that in the process 
all that may be left is a string of the four operators above. This proves the completeness of 
this set. 

Furthermore, in any string of S + , J + , A + , and A + , the order is unimportant up to the 
overall sign. This is because the commutators 

[S+,J + ] = [S+,A+} = iS+,A + } 
= [J + ,A + ] = [J+,i+] = [A+,i+]=0 (111) 

all vanish. A complete basis for the E = multiplet is thus formed by the states 

i w ,»i> = (s+r (jt (-t + r (*r io>. 
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These states are non- vanishing only if ha and are zero or one, as follows from the nilpotency 
of A + and A + , e.g. 



/ < J my m|J nj n| 

-E/IA/rL T C = -(^ + ) 2 = °- ( 112 ) 



Likewise, (A + ) 2 = 0. It is thus natural to break up the states into four ladders, 



(j+r io), ma) 

( S+ ) 



A*\e>, (ii4) 

'" 'j+)-i + |0), (115) 



(s + ) 

(S+)" F_1 (j + )"' i_1 ^/i + |0). (116) 
The allowed values of np and % in each ladder are essentially determined by the exclusion 



principle. In the first sector, Eq. (|113 ), one may add up to 2N Fermions to fill all N sites with 
both spin species, so that n F — 0, . . . , N. Since any number of Bosons may live at a site, the 
sum on n B = 0, . . . , oo is unbounded above. This result agrees with what we know from group 
theory: these states are simply direct product states of SU(2) and SU(1,1) representations. By 
computing 

S z \0) = \Y, (/it/»t - fnjli) |0> = -y |0>, (117) 

and recalling 5 _ |0) = 0, we recognize that |0) is the lowest weight state of a spin s = N/2 
SU(2) ladder, which has the well-known degeneracy 2s + 1 = N + 1. This agrees with the 
assignment S z (n F ) = -N/2 + n F . Likewise, J~\Q) = 0, J 2 |0) = f |0), and J z (n B ) = N/2 + n B . 
Thus |0) can also be regarded as the lowest weight state of a representation of SU(1,1), with 
spin s = N/2. This representation is irreducible, unitary and infinite-dimensional, where the 
third property is necessitated by the second one, since SU(1,1) is non-compact. 

Application of the Pauli principle determines the allowed values of n F in the other cases. 
In the second and third sectors, Eqs. (|114 - 115)) , one Fermion is already present, so at most 



2N — 1 may be added. Since they come pairwise, we must restrict n F — 0, . . . , N — 1 in these 
cases. Finally, in the fourth sector, Eq. ( |116|) , 2N — 2 Fermions may be added, suggesting that 
we need only n F < N . However, if n F = N, there are iV Fermionic pairs present, and the 
Fermionic content of the state is uniquely determined. That is, the full state must be a direct 
product of a Bosonic part with the fully occupied Fermionic state. This implies that 

S+) N (J+Y B |0> oc (S^ 1 (j+) nB - 1 A+i+|0>. 



We therefore choose n F — 1, . . . , N — 1 in the last sector. 
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5.3 Supersymmetric quadruplets 

To calculate traces over the E = multiplet, it is beneficial to organize the ladders of states in 
a slightly different way. To that end, consider the Fermionic operators 

Qa = £&L/na, (H») 
n 

which are particular combinations of the global currents i7tot> anc ^ belong to the set of twelve 
raising and lowering operators of which we have utilized only S + , J + , A + , and A + so far. The 
set {Qa,Q a ,nf a ,nb a } where a =| or a —I, closes under the graded commutator (i.e. the 
anticommutator for two Fermions, and the commutator in the other cases), and thus forms a 
subalgebra of the full algebra of global currents. This subalgebra is denoted H = u(l|l)©u(l|l). 
(The direct sum is over spin.) The significance of 7i derives from the fact that its generators 
commute not only with Hp but also with H v . Thus Hjj + H v is a Casimir invariant for 7i. 
Schur's lemma then states that Hp + H v is a multiple of unity on every irreducible multiplet 
of Ti. Let us therefore construct these multiplets. 

We interpret Qi, Qi as lowering operators, and Q-p Q± as raising operators for TL. A lowest 
weight state, |LW), for H obeys 

g T |LW) = Q i |LW) = 0. (120) 
To solve these equations for |LW), we take note of the commutation relations 

(121) 
(122) 
(123) 
(124) 

All other graded commutators between the two sets {S + , J + , A + , A + } and {Q^ : Qj., Qp Q^} 
vanish. Using the above commutation relations and making the ansatz 

ilw) = a(s+y(j+y>\ 0) 

+ e(s+y'- l (.ry"->A+A+\o), (125) 

one finds that the lowest weight condition can be satisfied with a = 1 and (5 = —n\. We denote 
the resulting lowest weight state by |n 1 n 2 ,00). 

Application of the raising operator to |n 1 n 2 , 00) produces 

Kn 2 ,10> = (^) ni (j+) n2_1 ^|0) 

cx Q^n^OO), (126) 

while application of gives 

1^^,01) = (5 + ) ni (j + ) n2_1 i + |o) 

oaQ l | n^, 00). (127) 



[Q h s + ] 


= A\ 


[Q P A + ] 






= A\ 


[Ql:A + ] 


= J\ 


[Q V J + ] 


= A\ 


[Qh A+ ] 


= J\ 


[Qi,s + ] 


= ~A\ 


[Qi,A + ] 


= -s + 
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When acting with the raising operator Q+ on \n x n 2 , 10), or with Q, on |%n 2 , 01), we get zero 

— 2 — 2 

from = = Q< due to the Fermionic nature of these operators. A non-vanishing state 
results from applying to 1%^, 10): 

1^,11) = (S+) ni+1 (J+Y 2 - 1 10) 

+ (n 2 -l) (Vf (j+)*"Vi+|0> 

oc QJn^, 10). (128) 

Exactly the same state is obtained by applying to \n l n 2 , 01). The state Eq. fll28[) is annihi- 
lated by the action of both and Q^, which again follows from Q 2 = = Q 2 . Thus \n x n 2 , H) 
is a highest weight state for H. 



t t t t t 




Figure 2: Illustration of the organization of states into 7i-quadruplets, shown here for the case 
N = 5. Each quadruplet contains a state with ri\ Fermion pairs and n 2 Boson pairs, plus 
three partners with equal "spin" but differing Boson/Fermion content. The n\ = 0, n 2 = 1 
quadruplet is shown as an example. The vacuum state, with n\ = ri2 = is a unique 7i-singlet. 

The states |n 1 ra 2 , [iv) (//, v — 0, 1) so constructed are seen to be in one-to-one correspondence 
with the previous set Eqs. ( |1 13 - 1 16 ) . What we have thus achieved is an organization of the 



zero energy sector into irreducible H-quadruplets (see Fig. 2), with the exception of only the 
vacuum |0), which does not fit into the above scheme but figures as a separate 7^-singlet. Each 
quadruplet consists of two Bosonic (fiu = 00, 11) and two Fermionic (fiu = 01, 10) states, 
reflecting the supersymmetry of the formalism. Recalling the earlier discussion based on the 
exclusion principle, we see that the allowed quantum numbers for the quadruplets are rii = 
0, 1, . . . , N — 1 and n 2 = 1, . . . , oo. The pair of quantum numbers n\ = n 2 = corresponds to 
the H-singlet vacuum state |0). 

By construction, the quadruplet states |n 1 n 2 , \iv) taken together with the vacuum, constitute 
an orthogonal basis of the zero energy sector. As it stands, they are not normalized to unity. 
While it is not difficult to include the correct normalization factors, there is no need to do that 
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here. For our purposes, all we require is that we be able to calculate traces, and for that linear 
independence of the states is entirely sufficient. 



5.4 DOS Correlation Function 

We are finally in a position to calculate correlation functions. As a check, consider first the 
partition function within the zero energy multiplet. Since the full partition function equals 
unity for all values of the parameter d and approaches the zero mode contribution as d — > oo, 
the zero mode contribution must itself be exactly normalized to unity. And, indeed, 

Z = Tr (-l) NF e- pH = 1, (129) 

where the subscript signifies the restriction to the zero energy sector. The 1 on the right- 
hand side stems from the vacuum. All other contributions cancel as H is constant on each 
7^-quadruplet, and the two Fermionic states of a quadruplet come with a minus sign relative 
to the two Bosonic ones, cf. Sec. [3i)i Since H v , by Eq. (|2~4"|), just counts the total number of b 
Bosons and / Fermions, its action on the quadruplet states is 

H^\n^n 2 ,iw) — 1r\{ri\ + n 2 )|n 1 n 2 , (130) 

Now consider the DOS correlator, Eq. ( |101| ). Defining 

G± n {oj) = G±(n, n; x, x; uj) (131) 

and recalling Eq. flTT|), we get 

p(w/2) = /t£ [G-»(w/2) - G+niu/2)). (132) 
2m n 

The product therefore becomes 

4n 2 p(u;/2)p(-Lj/2)/f3 2 = 



E 



G +n {u/2)G- n ,(-u/2) + G- n (u/2)G +n ,(-u/2) 



-G +n (u/2)G +n ,(-u/2) - G^ n (u/2)G- n ,{-u/2) 



(133) 



To proceed, we perform an ensemble average, which may be split into connected (cumulant) 
and disconnected pieces. To calculate the disconnected terms, we use the symmetric operator 
representation 

G+n = -\{fnA-fy^) =</It/»T- < A ( 134 ) 
G-n = -\(fljnl- fulfil) =-^1+^ (135) 



Therefore: 



G +n {u/2)G^{-u/2) + G_ n {u/2)G +n ,{-uj/2) 
-G +n {u/2)G +n ,{-uj/2) - G_ n (uj/2)G_ nl (-uj/2) = 1, (136) 
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where G± n = [G± n ]ens = These disconnected terms thus cancel the — p 2 term in the 

definition ofC(w), Eq. (pip . Then 



n,n — 1 



2 AT 

, G +n (u;/2)G_ n ,(-a;/2) + G^u/2)G^-u/2) 



-C +n (a;/2)G +n ,(-a;/2) - C_ n (a;/2)G_^(-a;/2) 



(137) 



where the subscript c indicates the connected average. 

It is well known that the last two terms do not contribute to the correlator .0 Using this 
simplification, we arrive at the form 



? 2 v 
2^ 



C(u) = ^Re £ [G +n (^/2)G_ n ,(-^/2)] cnSiC . (138) 

n,n'=l 



As with the diffuson, this may be calculated in several ways, equivalent by supersymmetry. 
We arbitrarily choose the Fermion representation 



2tt 
2tt 2 



Re(n n n n ). (139) 



The quadruplet states |ra 1 ra 2 , pv) are eigenstates of n^Ufi with eigenvalues n 2 , n 1 (n 1 + 1), 
(ni + l)nx, and (ni + l) 2 for /xt/ = 00, 01, 10, and 11 respectively. Therefore the trace of n^Ufi 
over the quadruplet with quantum numbers n\ , n<i is 

n\ - 2m(ni + 1) + (ni + l) 2 = 1. 

Using this to evaluate the expectation value Eq. (|139|) as a trace against the Boltzmann weight, 
we obtain 

r>2 oo N—l 

C(u) = -^Re E e - /3(2 " + ^ ){n2+ni) . (140) 
2 71 " „ 2 =i„ 1= o 

While this sum may, of course, be evaluated exactly, it is simpler to take the continuum limit, 
valid for f3u <C 1. Measuring energy in units of the level spacing A = 2tt/N/3, we define 
uj = Nf3u /27r, and rescaled integration variables x = n 2 /N, y = n\/N. Equation Eq. ( |140| ) 
then becomes 

C(u) = — ^-Re / dx dy e ~ 2m ^ x+y \ (141) 
2% 2 Jo Jo 

After division by p 2 = A~ 2 , and omission of a 5(ui) term, this result takes the universal form 

CM/P 2 = (142) 

which is well-known from random matrix theory [Q. Similar manipulations lead to identical 
results when C(o;) is evaluated using the purely Bosonic or mixed expressions analogous to 
Eq- Q. 

2 If they are calculated directly using a modified supersymmetric functional, the resulting Hamiltonian has a 
unique ground state even for r\ = 0, indeed leading to this conclusion. 
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5.5 Extension to parametric correlations 

The above calculation can be extended to the case where the random potential V n (x) depends 
on a parameter, A: 

= E / dxV n (x,\)^ n , (143) 

n J 

V n (x, A) = V n {x)-XW n (x). (144) 

The extra term XW n (x) will be called the "background potential". For simplicity we take 
it to be independent of x: W n (x) = h n . (It can be shown that the final result ( |158| ) remains 
essentially unchanged when this assumption is lifted.) We assume that there are no correlations 
between the background potential Xh n and the random potential V n (x). For the moment the 
other statistical properties of the numbers h n are left unspecified. 

The object of interest is the parametric density of states correlation function 

C(^A) = [p(^A)p(0,0)] cns -p 2 . (145) 

By the same steps as in the previous section, its calculation reduces to that of 

E [G+n(w/2, A/2)G_ n ,(-u,/2, -A/2)] enSiC . 

n,n' 

The calculation bears much similarity to the one we did for A = 0, and we will now focus on 
the new features arising for A ^ 0. 

The background potential gives rise to an additional piece in the Lagrangian, 

£ n -> C n + l -Xh n (0*cr 2 n + $ n Vzipn) , (146) 
and in the spin Hamiltonian, 

H^H + XH M , H M = iY,K{S z n + J z n ) ) (147) 

n 

which acts like the coupling to an imaginary magnetic field. Note that adding XHm is the same 
as substituting uj — > uj + Xh n . To prevent the background potential from causing an overall 
shift of the frequency, we shall require that J2n h n vanishes identically for every realization of 
the disorder. 

What is the effect of H M inside the zero energy multiplet? By the requirement J2 n hn = 0, 
the zero mode approximation S* + J% —> N' 1 J2n(^n + J*) for H M gives a vanishing result to 
linear order in A. To see the effect of the background potential, we must go to higher order and 
do a calculation similar in spirit to that performed by Simons and Altshuler[p3| for the case of 
time-reversal invariant disordered metallic grains. What the perturbation Hm does is to couple 
the zero energy multiplet with the one-magnon sector. By second-order perturbation theory, 
this coupling results in an effective Hamiltonian, H\, acting on the states of the zero energy 
multiplet: 



H x = Hu\-^2^}H M , (148) 
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where Hk is the projector onto the one-magnon sector with momentum k. In the appendix it 
is shown that for N large, H\ reduces to 



H, = ^~, 7=£e^, (149) 

where 

C+- = S + S- + J + .r + A + A~ - A + A~, (150) 

and 

^ = ^E e ^n (151) 

n 

is the Fourier transform of the background potential. In order for second-order perturbation 
theory to apply, the change in energy due to the perturbation must be much smaller than the 
smallest energy denominator. Because the gap for single magnon excitations is of the order of 
D/N 2 , and H\ will be seen below to shift the energy by an amount of the order of jN 2 , the 
condition 

7 < D/N A (152) 

on 7 is sufficient for the perturbative formula (|148|) to be valid. To give this result some added 
physical meaning, consider the simple case in which the numbers h n = h n — h, where the h n are 
chosen to be independent, identically distributed random variables with zero mean and unit 
variance, and h = N^ 1 Y^ n hn is chosen to maintain zero spatial average. The ensemble average 
of the quantity \hk\ 2 is 1/N (for k ^ 0), and 

Wens= N 2 ^'Dk 2 = 12D (153) 
in the large N limit. More generally, we are led to define a dimensionless, order one, 70 by 

A 2 

7 == 7j7o- (154) 

In the particular example here, comparison of Eqs. ( |153| - |154D gives [7o] e ns. — 1/12. The condition 
on the smallness of 7 then becomes 

A < D/N 2 , (155) 

i.e. the typical strength of fluctuation of the background potential must be much smaller than 
the one-magnon gap. 

We are now ready to compute the parametric density of states correlator, C(u,X). From 
Eqs. (JL21[ - |124|) and similar commutation relations involving the lowering operators S~ , J~ , A~ 



and A~ , one finds that C + ~ commutes with the Fermionic charges Q a and Q a , which means 
that C + ~ is a Casimir invariant for TC = u(l|l) © u(l|l). Therefore C + ~ must be proportional 
to unity on the quadruplets: 

C +_ |n 1 n 2 , jj,u) = T(ni, n2) |n 1 n 2 , /iz/). (156) 
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To calculate the eigenvalue r(ni,n2), we may work out the action of C + on any one of the 
four states (/i, v = 0, 1), say |ri 1 n 2 , 01). A straightforward calculation gives 

r(ni,n 2 ) = n 2 (N + n 2 -l) 

+ m(JV-m-i). (157) 

By including the contribution from H\ = / -fC~ i into the Boltzmann weight, we obtain for the 
parametric density of states correlator, 



C(u, A) = -^Re J2 E e~ i/Mn2+ni) 



2 oo N-l 

n 2 =lni=0 

xe -/3A 2 D- 1 7 or(n 1 ,n 2 )_ /-^ 



For (3N\ 2 D~ lr )o <C 1 or, equivalently, A 2 /} _1 7o <C A, we may again take the continuum limit 
and replace the sum by an integral. By performing a simple but revealing substitution of 
integration variables, we get 

°(^ x ) = f I e <o(p-,)-AV-« 2 ) 4p^? n 59 ) 

P 1 J \p\>k J \q\<-* (2tt) 2 ' 
with the rescaled parameter A given by 

N 2 

A 2 = 0— X 2 lo/(2ir) 2 . (160) 

Upon making the identification of to with position and A 2 with imaginary time, the expression 
for C(u,X) is seen to coincide with the dynamical density- density correlation function of a free 
Fermi gas. For the case of small metallic particles, this coincidence was recently pointed out 



by Simons et al . [ 2 1 



6 Conclusions 

6.1 Summary of results 

We have studied the zero-temperature properties of the surface states of a three-dimensional 
layered integer quantum Hall sample. Previous work|)|, |10| showed that for weak disorder in 



the thermodynamic limit this system is a dirty chiral metal, with finite conductivity along the 
field axis at zero temperature, and ballistic transport transverse to it. 

In this paper, several new results are obtained and placed in a more powerful concep- 
tual framework. By constructing an appropriate generating functional for disorder averaged 
correlations, we showed that this system can in fact be viewed as a kind of supersymmetric 
(SUSY) ferromagnetic spin chain. Like for an ordinary ferromagnet, the ground state and single 
magnon excitations can be obtained exactly. This allows an exact calculation of the diffuson (or 
density-density correlator). The conventional z = 2 dynamical scaling for ferromagnetic spin 
waves corresponds to the combination of finite conduction (transverse diffusion) and ballistic 
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transport discussed above. In the thermodynamic limit, therefore, we have demonstrated the 
absence of localization in the 2d chiral metal, even in the limit of strong disorder. 

In addition, the quantum formulation of the SUSY ferromagnet provides a novel method to 
study the mesoscopics of the chiral metal. As an example, we have determined the correlation 
function of the density of states at two energies and for two slightly different realizations of 
disorder, the so-called parametric correlations. As expected, in the ergodic regime this takes 
on the universal form appropriate for the unitary ensemble, i.e. a free Fermion density-density 
correlation function . 



6.2 Connection to the non-linear a model 



It is useful to connect the approach taken here to the more conventional NLcxM method pA 
This connection may be made quite explicit using the SUSY coherent states developed in the 
appendix. We will construct a coherent state path integral representation for the partition 
function Z. We begin by considering the problem of decoupled edges, i.e. with D = 0. Then 

Z = T[Zn, (161) 

n 

where 

Z n = STre-P Hn . (162) 

The single super spin Hamiltonian for A = is determined by the condition H v = J2 n H n , so 
by Eq. ©, 

H n = 2 V (S z n + J z n ) . (163) 

The coherent state path integral can then be simply obtained using the resolution of unity, 
Eq. Q174]) , and the formula for the supertrace, Eq. ( |175| ). One finds 



Z n = / D(Z n ,Z n )e- s "^\ (164) 

where the single mode action is 

rP rp 
S n [Z n ] = dr (Z n \d T + H n \Z n ) = / drC n . 
Jo Jo 

The single spin Lagrangian is given by 

£n = (Z n \d T + H\Z n ) 

= -STr(l + Z n Z n y 1 Z n (d T + 2 V )Z n . (165) 

This expression can be put into a simpler form by defining the 4x4 supermatrices 

*=(-k t). A -(J -J < 166 > 

and 

Qn = g n Ag- 1 - (167) 
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In terms of these fields, the single spin Lagrangian becomes 



C n = isTr (Ag-%g n ) + |STr (AQ n ) . (168) 

The first (dynamical phase) term is of the Wess-Zumino type, and cannot be written in a 
globally non-singular form in terms of the Q field. A Wess-Zumino term also occurs in the 
coherent state path integral for an ordinary SU(2) spin|23[], and could have been expected here 
on general grounds. It is necessary to obtain the z = 2 dynamics appropriate for a ferromagnet. 
Inclusion of such a term allows us here to obtain a true NLaM formulation, in contrast to the 



earlier Q- matrix formulation in Ref.|10[, in which only an expansion in the ordered (metallic) 
phase was determined. 

Eq. ( |168| ) holds for a single super spin. It is now straightforward to include the exchange 
coupling H D to obtain 

Z = J D(Z,Z)e~ s (169) 

where 

S= f drJ2(C n + C^ n ), (170) 
Jo 

and 

£int,„ = jSTr {Q n Qn + l) ■ (171) 

Eq. ( |17Up is the full NLaM action for the 2d system. The universal level correlations, usually 
discussed in this context, are properties of the zero-dimensional quantum limit, in which it 
is appropriate to make a "zero mode" approximation, neglecting spatial and temporal (i.e. 
t = x) variations of Q. The zero mode theory then becomes a single super-integral over Qq. 
This formulation can be connected back to the operator one of the text, by noting that the 
resulting integral is simply another representation of the restricted trace over the zero energy 
multiplet of completely polarized super spin states. 

6.3 Questions and open problems 

We conclude with a discussion of various questions which remain open to future investigation. 
One interesting issue is to understand in more detail the nature of the crossovers between the 
three mesoscopic regimes described in Section V. In the language of the ferromagnetic super 
spin chain, which has finite length N and is at finite temperatures {(3 < oo), the three regimes 
correspond to: (i) a "zero- dimensional" regime in which all (finite wavevector) spin-waves 
are absent, being too costly in energy, (ii) a "Id diffusive" regime in which the spins behave 
classically but are still ferromagnetically ordered, and (iii) a "localized" regime in which the 
spin chain length exceeds the ferromagnetic spin correlation length. For ordinary (non-super) 
ferromagnetic spin chains, recent progress has been made in computing various scaling functions 



connecting these regimes |27]]. The apparent generality of the technique suggests that progress 
might likewise be made for the super spin chain. |]29[ 

Also of interest in terms of crossover behavior are the scaling properties of the wavefunctions 
in various limits. Although we do not expect interesting scaling properties in the thermody- 
namic limit (because the 2d system is truly a stable (chiral) metal), interesting possibilities 
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arise in the intermediate mesoscopic regime. Specifically, in the "Id diffusive" regime, with 
L = y/D(3 < iV C ( = DP, the electron motion is ergodic on frequency scales uj ~ A. 
Nevertheless, there is a breakdown of the zero-mode approximation, which requires Lq ^> N. 
We then expect the wavefunctions to become very broadly distributed random variables]^. It 



would be interesting to see if the moments exhibit multifractal scaling in this regime, as well 
as to look at scaling of multi-point correlators. 



Thirdly, it would be interesting to calculate the conductance fluctuations |30|| , which would 
involve applying appropriate boundary conditions to the ends of the super spin chain and cal- 
culating an eight-Fermion (or eight-Boson, etc.) correlator (product of four Green's functions). 
Preliminary work using diagrammatic techniques by Mathurj3l| suggests rather interesting 
crossover phenomena in the variance of the conductance. 

Lastly, these methods may be useful to connect with earlier work on chiral classical wave 



propagation in a disordered medium. As pointed out in Ref. P2fl , the Schrodinger equation 
for electrons in the surface sheath may be viewed as a classical wave equation for "directed" 
propagation. It is then of interest to study the spreading and deflection of a point source at an 
initial x = 0,n = to some larger x. The beam width is closely related to our diffuson, but 
(like the conductance fluctuations) the motion of the beam center involves the product of four 



Green's functions, but may also be tractable [|3T . 
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Appendix 

This appendix is in two parts. In the first part we are going to develop some of the mathematical 
structures underlying the coherent state path integral for the super spin Hamiltonian. In the 
second part we will derive the expression (|149|) for the effective zero mode Hamiltonian H\, 
which results from treating the background potential by second-order perturbation theory. 

Our first item will be to discuss the group of canonical transformations of the Bose and 
Fermi operators. For that purpose we introduce a more economical notation, by setting 

cb T = B h c F i = F h cbi = E>i, c F [ = F;, 
c B \ = B h c FT = F\, cbi = B b c Fl = F\. 

Now let A, B, C and D be complex 2x2 supermatrices, i.e. 

A=( A / F A / B ) etc., 
\A BF A BB J 

where App and Abb are complex numbers, while A F b and Abf are Grassmann numbers. 
Consider then the transformation 

r n r n A -I- r n C 
C al ^ c b1 n ba T C b ^U ba , 
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where the superscript n = 1, N indexes the chiral modes and, here and below, the summation 
convention is used. By setting {c a } = {cf^, csj, cfj., Cbj,}, we write this in the abbreviated form 

A B 



'a 



If g has an inverse, we can also transform the annihilators, by 



c -> (g ' ),. 



When rearranging products, we follow the convention that Grassmann numbers not only an- 
ticommute among themselves, but also anticommute with the Fermionic operators. It is then 
easy to see that the transformation c™ — > (g~ 1 ) a/3 c /3 , c™ — > c^g^ a , is canonical, i.e. preserves the 
graded commutator. 

Canonical transformations have an inverse, and a succession of two canonical transforma- 
tions is again canonical. Therefore such transformations p3|| form a group, which in the present 
context is the Lie supergroup Gl(2|2), obtained by complexifying U(2|l, 1). The elements g of 
Gl(2|2) can be written in the form of a Gauss decomposition, 



9 



C DJ ~ \0 l) \ hJ\Z 1 



where, from an easy calculation, 

Z = BD-\ hi = D, 

Z = D^C, hf = A- BD^C. 

Every canonical transformation g can be realized by a Fock space operator T g that is ob- 
tained by exponentiating some bilinear cc, and acts by c™ — > T g c 7 ^T~ 1 , c™ — ► T^c^T" 1 . In detail 
the correspondences are as follows: 



1 



T /l 0\ = exp (c^Z^) , 
\Z 1 



T fe A 0\ = ^Pfe^ 4 



1 



T fl \ = eX P{ C al D abC H 
,0 



By construction, the correspondence g —>■ T g defines a representation of Gl(2|2) on Fock space. 
Therefore every relation that is valid for supermatrices g, and uses no more than the Lie 
supergroup structure of Gl(2|2), also holds for the Fock space operators T g . In particular, by 
applying the Gauss decomposition to the product 

1 0W1 Z 
Z 1 J VO 1 
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we get 

= exp(c^[Z(l + ZZ)- l ] ab cl) 

x exp (c^[ln(l + ZZ)} ah c n hi - c^[ln(l + ZZ)} ah c n b . 



e W (c n al [Z(l + ZZ)- l [ 



ab C b*\ 



which is called a "disentangling" identity, since it moves the Gl(2|2) lowering operators all the 
way to the right, and the raising operators all the way to the left. 



Consider now the generalized coherent states p4 



\Z) = exp faZtty |0)SDet(l + ZZ)^ 2 
(Z\ = SDet(l + ZZ) N / 2 (0\exp(clZ ba c n a 

The superdeterminant of a 2 x 2 supermatrix, 



SDet 



M FF M FB 
M BF M bb 
= {M bb - M BF M FF - l M FB ) /M FF 

= M BB j (M FF - MfbMbb^Mbf) , 

satisfies InSDetM = STr InM, where STrM = —M FF + M BB =: Ea(^) lal M aa , with \a\ = if 
a = B (Boson) and \a\ = 1 if a = F (Fermion), is the supertrace. Using this, the disentangling 
identity, and the relations 

SMfl°> = °> c n aiCbi\0) = 5 a b(-l) lal+1 N\0), 

it is easy to check that the generalized coherent states are normalized, 

(Z\Z) = 1. 

Another useful way of writing the generalized coherent states is 

\Z) = T g \0), (Z\ = (0\T-\ 



where 



and 



1 Z\ ((l + ZZ) +l ' 2 \ / 1 

lj\ {l + ZZ)- 1 ' 2 ) \-Z 1 

(i + zz)- 1 ' 2 zii + zzy 1 / 2 
-z^i + zzy 1 ' 2 (i + zz)- 1 ' 2 



{i + zzy 1 ' 2 -zii + zzy 1 ' 2 
z(i + zz)- 1 / 2 (i + zz)- 1 / 2 



(172) 



(173) 
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The zero energy multiplet of Sec. coincides with the space of states obtained by acting 
repeatedly with the Gl(2|2) raising operators c£+c£r on the vacuum |0). In|]35| it was shown 



that, if D(Z, Z) = Dg H is the uniform superintegration measure on the supermanifold G/H = 
G1(2|2)/G1(1|1) x Gl(l|l) parameterized by the coherent states \Z) = T g \0), the unit operator 
on the zero energy multiplet can be resolved by 

1 = J Dg H T g \0)(0\T- 1 = J D(Z,Z)\Z)(Z\. (174) 

(Note that the conventions of that reference differ from ours by Z — > —Z and the exchange of 
the Fermionic and Bosonic sectors.) The integral is over 

Z FF = Z FF} Z BB = —Z BB) 

with the integration domain being 

< \Z FF \ 2 < oo, < \Z BB \ 2 < 1. 



These relations mean[T8|] that the variable Z FF is a complex stereographic coordinate for the 
two-sphere S 2 (as is well-known for SU(2) coherent states) while Z BB parameterizes a two- 
hyperboloid H 2 (as is appropriate for the non-compact group SU(1,1)). Of course Z BF , Z FB , 
Z BF , and Z FB are Grassmann variables. 

Given the resolution of unity ( [174] ), we can convert supertraces over the zero energy multiplet 
into integrals over the generalized coherent states: 

(O) = Tr (-l)^O = J D{Z,Z)(Z\0\Z). (175) 

The coherent state path integral for the partition function of a super spin system with Hamil- 
tonian H, is now obtained in the usual manner by inserting resolutions of unity between in- 
finitesimal imaginary time slices of the Boltzmann weight exp(—/3H). The Lagrangian of the 
resulting path integral is £ = (Z\d T + H\Z), whose explicit form is easily calculated by using 
the definition of the coherent states and the disentangling identity. For example, 

(Z\d T \Z) - d T lnSDet{l + ZZ) N/2 

= SDet(l + ZZ) N (0\ exp(c^ 6a < T )«9 T exp^^c^O) 

d _ 
= SDet(l + zz ) N g^ s - SBet ^ + zz + sZd T Z)- N 

= — iVSTr [(1 + ZZ)- 1 Zd T Z . 

The Wess-Zumino (or linear in d T ) term in ( |165| ) is obtained on setting the number of chiral 
modes TV equal to one and omitting a total r derivative. 

Below we will need an explicit expression for the coherent state expectation value of the 
operator C + " appearing in the formula (|149[) of H\. To calculate it we write 

C+- = S + S- + J + .r + A+A- - A + A~ 
= S~S + + J-J + + A-A + - A~A + . 
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The second equality is invalid for both SU(2) and SU(1, 1), since S + S ^ S S + and J + J ^ 
J~J + , but it does hold in the present result of cancellations due to supersymmetry. 

Using the second form of C + ~ we easily find 

SDet(l + ZZ)~ N (Z\C + -\Z) 

= E(-l) |d|+1 (0|exp(^^ a c: T )exp(^ T Z a ^)|0) 
aZj cd oZj dc 

= E(-l) |6|+1 -^SDet(l + ZZ)-", 



which yields 



(Z\C + -\Z) = -N 2 STr [ZZ(l + ZZ)- 2 \ + 0(N l ). (176) 
We now turn to the effective zero mode Hamiltonian 

and show that it reduces to ( |149| ) in the large- N limit. The first step is to construct the 
projector on the space of momentum-^ single magnon excitations, life. Magnons can be created 
on top of the vacuum or any other state of the zero energy multiplet. We define the single 
magnon creation and annihilation operators by 



^ab(k) — E e ^ C al C bp 



\ab,k) = -j=J+(k)\0). 



and consider the states 



The corresponding projection operator is 

7r fc = ^|a6 ) A;)(a6,fc| = l^(-l)l b l + V a |W|0)(0|J 6 ;(A;). 

ab ab 

It is easy to verify that, with the normalization and sign factors chosen, tt^ satisfies — i?ki 
and is a singlet with respect to the (complexified) algebra TC = gl ( 1 1 1 ) © gl(l|l) generated by 
the zero mode operators Q^, Qi, Qp Q,, n^, n/;, n^, and n^. The states \ab,k) are single 
magnon excitations built on the vacuum. Other one-magnon states (not based on the vacuum) 
are generated by applying a global rotation, 

\ab, k) — > T g \ab, k). 

The projector 7T& transforms as 

The projector on the entire one-magnon sector is obtained by averaging over all global rotations, 

U h = J Dg H TgTTf.Tg" 1 . (177) 



31 



The integrand is well-defined as a function on the coset space G/H since, by the 7i-singlet 
property of 7r fc , 

TghFkTgh = Tg^h^k^h T g = T g 7T k T g 

for h G H = Gl(l|l) x Gl(ljl). 

To go further, we need to take advantage of the large- A?" limit. Let us evaluate the expec- 
tation value of H M Tl k H M in a coherent state \Z) = T g \0), 

E:=j Dg' H (O^H^ir^H^O). 

For large N we expect the dominant contributions to the g' integral to come from the imme- 
diate vicinity of g — g(Z, Z), so we substitute variables g' = gexpX (with X subject to the 
transversality condition AX + XA = 0) and approximate Dg' H by its linearization DX at g. 
The integration over X is trivial since 

(OlT-XxpxlO) = (0|T cxpX |0) ~ e -^ STrf2 /4 _ S{x) 

reduces to a 5-function in the large- N limit. Thus, doing the g' integral to leading order in the 
large parameter N, we get 

E = (0\T- l H M T g n k T g - 1 H M T g \0). 

The large- X limit is semiclassical in nature, and it is not hard to see that all off-diagonal matrix 
elements of H\ become negligible for N — > oo. Therefore, it is in fact sufficient to consider the 
diagonal ones, which is what we are doing. 

The next step is to compute the rotated perturbation T g -1 iJ M T g from the expression H M = 
i En,a Mc" T c" T - c^)/2. For this we set 

(A B\ _! (A B\ 

9= {c d)> 9 ={c 

and use the transformation laws (see the beginning of this appendix) 

rp-l-n rp _ - n a i_r n r< 
1 g L a\ 1 g ~ L a'\ IX a'a + C a'l U a'ai 

T g lc a\T g = A aa ,c a a ,^ + -B aa /C™^, 

which yield 

T 9 - l H M T g = \- ^2 h n (c™(AB — BD)c™ 

" n 

+c n l (CA- DC)tf + ...). 

On taking the matrix element with the one-magnon creation operator, we obtain 

jN 

m; l H M T g J+{k)\V) = —h k {-lt m+ \CA-DC) ba . 
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Similarly, 

(OUbaWT-'HuTslQ) = —h- k (-l)^ +1 (AB - BD) ab . 
We now switch notation from g to Z, Z. From ( 172 , 173 ) we read off the relations 



CA-DC = 2Z(l + ZZ)-\ 
AB-BD = (l + ZZ)- l 2Z, 

whence the coherent state expectation value of H\ takes the form 



(z\h x \z) = -nJ2 x STr zz ^ 1 + ZZ Y 



Comparison with (|176[) gives H\ = 7C + , as claimed in ( |149| ). 
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